Condensation-Driven Aggregation in One Dimension 
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We propose a model for aggregation where particles are continuously growing by heterogeneous 
condensation in one dimension and solve it exactly. We show that the particle size spectra exhibit 
transition to dynamic scaling c(x,t) ~ i - ' 5 <)>(x /t z ). The exponents 3 and z satisfy a generalized 
scaling relation 8 = (1 + q)z where the value of q is fixed by a non-trivial conservation law. We have 
shown that the value of (1 + q) is always less than the value 2 of aggregation without condensation. 

PACS numbers: 61.43.Hv, 64.60.Ht, 68.03.Fg, 82.70Dd 



I. INTRODUCTION 

The formation of clusters by aggregation of particles 
is a characteristic feature of many seemingly different 
processes in physics, chemistry, biology and engineer- 
ing. Examples include aggregation of colloidal or aerosol 
particles suspended in liquid or gas [l], [1, 9) polymer- 
ization antigen-antibody aggregation y| , and cluster 
formation in galaxy Q. This wide variety of applica- 
tions has resulted in numerous studies which reveal that, 
when chemically identical particles aggregate, almost al- 
ways scale invariant clusters emerge. Note that due to 
the non-equilibrium nature of the aggregation process 
the standard theoretical framework developed for equi- 
librium statistical physics is found redundant. However, 
the application of stochastic theory is found to be increas- 
ingly useful in capturing a wide class of non-equilibrium 
phenomena. 

Typically, the non-equilibrium systems are described 
by the rate equation approach, having the form of a 
master equation which is often governed by some con- 
servation principle. The Smoluchowski equation for the 
kinetics of irreversible aggregation is one such example, 
where the distribution function c(x, t) of particle of size 
x at time t evolves according to the following integro- 
differential equation 0, [1] : 

dc(x,t) . . [°° r „ w . , 

— — — = -c{x,t) J K{x,y)c(y,t)dy 
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dyK(y,x-y)c(y,t)c(x-y,t). (1) 



Here, the kernel K(x,y) is symmetric with respect to its 
argument and it determines the collision time in which 
a particle of size x collides with another particle of any 
size y and they merge into a particle (aggregate) of size 
(x + y). The first term on the right hand side of Eq. ((T|) 
describes the loss of a particle of size x due to merging 
of particles of size x with another particle of any size, 
while the second term describes gain of x due to merging 
of particles of size (x — y) with y. The Smoluchowski 
equation has been studied extensively for a large class of 
kernels satisfying K{bx, by) = b x K(x, y). The homogene- 
ity exponent A is shown to play a crucial role in classify- 



ing gelling and non-gelling models. For instance, A < 1 
describes the non-gelling model whose dynamics is gov- 
erned by the conservation of mass principle, and A > 1 
is the gelling model that describes the gelation transi- 
tion accompanied by the violation of mass conservation 
law d,[T3|. Note, however, that, despite the seemingly 
simple structure of Eq. (Q]) , it is solved exactly for non- 
gelling model, only for a constant kernel A = 0, and that 
solution was given more than one hundred years ago by 
Smoluchowski himself. Finding another exact analytical 
solution of Eq. (fTJ) for the non-gelling model (A < 1) 
therefore still remains an open challenge. 

One of the reason why the Smoluchowski equation was 
so successful is that it has provided much of our theoreti- 
cal understanding about dynamic scaling and associated 
exponents which are in aggreement, at least qualitatively, 
with those extracted from real experiments and numeri- 
cal simulations pT| . For instance, for A < 1 it was shown 
that the distribution function c(x, t) exhibits dynamic 
scaling, 



c(x, t) ~ t~ p (f){x/t z ) ; with z > 0, 



(2) 



in the long time (t — * oo), large size (x — > oo) limit where 
4>(0 is a scaling function whose argument £ = x/t z is a 
dimensionless quantity. The exponents (3 and z satisfy 
a scaling relation (3 — Oz with 6 — 2 which follows from 
the conservation of mass principle. A scaling form like 
Eq. @ is shared by an extraordinarily diverse range 
of other phenomena, not just in aggregation, e.g., sys- 
tems exhibiting self-organized criticality, cluster growth 
in driven diffusive systems, fragmentation processes, etc. 

. The ubiquity of this scaling form suggests the ex- 
istence of a common underlying mechanism which makes 
such seemingly disparate systems behave in a remarkably 
similar fashion. 

In addition to growth by aggregation, there exists a 
host of other mechanisms (e.g., condensation, deposi- 
tion, and accretion) whereby par ticles can grow contin- 
uously between aggregations [l4 flEI . ITil QjJ. For in- 
stance, aerosol or colloidal particles are often not stable 
but evolve via aggregation and condensation, leading to 
gas-to-particle conversion. However, when the concentra- 
tion of particles present is high and the super-saturation 
is low, the condensation is heterogeneous in nature since, 
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FIG. 1: Schematic representation of the model for a = 1. 



in this case, condensation takes place onl y on the existing 
particles without forming new nuclei [ull6j|. Otherwise, 
the system may have sufficient number of impurities, such 
as dirt or mist particles, which usually serves as potential 
nucleation sites on which condensation takes place, and 
the resulting process is known as homogeneous conden- 
sation. In the latter case, the gas starts condensing on 
such nuclei and thereafter it is counted as a new particle, 
which then may take part on aggregation processes with 
other particle in the system. However, in the present 
work we will consider only the former case - heteroge- 
neous condensation. Such condensation-driven aggrega- 
tion does play an important role in the formation of the 
size spectra that ultimately determine the various phys- 
ical properties of the aggregates. Motivated by this, we 
here propose a simple model which is defined in the next 
section. 

The rest of the paper is organized as follows. In Sec. 
II, we present the definition of our model including its 
algorithm. Interesting results from numerical simulation 
based on the algorithm are presented in Sec. III. In sec- 
tion IV, we propose a generalized version of the Smolu- 
chowski equation and its appropriate parameters to de- 
scribe the evolution of the distribution function c(x, t) 
following the rules of the model. In Sec. V we give an 
explicit time dependent solution for c(x, t) and in Sec. VI 
we give its scaling description to obtain various scaling 
exponents. Finally, in Sec. VII, we discuss and summa- 
rize the work. 



II. MODEL 

We assume that initially the system consists of a large 
number of equal sized chemically identical particles which 
are assumed to be immersed in the gas phase. These 
particles, while in Brownian motion, are continuously 
growing by heterogeneous condensation, leading to gas- 
to-particle conversion, and merge irreversibly with other 
similarly growing particles upon encounter. To further 
specify our model, we assume that the amount of net 
growth by condensation of a given particle between col- 
lisions, in the most generic case, is directly proportional 
to its own size. The algorithm for one time unit of the 



FIG. 2: ln[L] vs ln[t] (where t = 1/iV since the number of 
particles present ultimately determines how fast or slow the 
aggregation process should proceed) are drawn using simu- 
lation data from one realization. The solid lines represent 
theoretical predictions with gradient equal to 2a. 



model can then be defined as follows: 

(i) Two particles are picked randomly from the system 
which mimics random collision vis-a-vis Brownian 
motion in one dimension. 

(ii) The sizes of the two particles are increased by a 
fraction a of their respective sizes. 

(iii) Their sizes are combined to form one particle. 

(iv) The steps (i) — (iii) are repeated ad infinitum. 

While the model is rather simple to define, the results it 
offers are far from simple. In order to illustrate how the 
rules of the model work for mono-disperse initial con- 
ditions, we give a simple example in Fig. 1 using an 
evolutionary-tree-based approach. The state of the sys- 
tem at a given time is fully described by the numbers in 
the corresponding box. The evolution of these numbers is 
then described by a set of such boxes in successive times 
along the possible trajectory, e.g., 1 — > 3 — ^6 — > and so 
on. The first salient feature of this model is that the sum 
of all the numbers L in the successive boxes keep increas- 
ing continuously with t revealing that the conservation of 
mass law is clearly violated. Secondly, the numbers in the 
different boxes represent the size or mass of the aggregate 
and hence a given box can be characterized by a distri- 
bution function c(x, t) of particle of linear size x at time 
t. 



III. NUMERICAL SIMULATION 

In order to extract a couple of basic features of the 
model, we have performed extensive numerical simula- 
tions based on the algorithm (i)-(iv). Perhaps it is worth- 
while to recall the work of Falk and Thomas, who in 
1974 obtained the molecular-size distribution by simulat- 
ing the discrete version of Eq. ([5]) [HI . The most crucial 
aspect of the stochastic simulation is to find a way to 
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FIG. 3: ln[s] vs hx[t] are drawn using data collected from nu- 
merical simulation. Data points (denoted by symbols) for 
a given curve represent one realization of our simulation. 
The solid lines represent theoretical predictions with gradi- 
ent equal to 1 + 2a. 



define the time t. Note that the time necessary for a par- 
ticle to come into contact with another particle should 
depend on the number of particles present in the system. 
Indeed, ultimately it is the number density of particles 
present in the system which should determine how fast 
or how slowly the process should proceed, provided the 
collisions for aggregation of particles are independent of 
their size, and hence we define time t = 1/N. In Fig[2]we 
therefore have plotted ln[L(i)] against ln[t] for different 
a and found a straight line with slope equals to 2a. This 
implies that 



(3) 



which immediately confirms that, for systems describ- 
ing the condensation-driven aggregation, conservation of 
mass principle is always violated. In Fig. [3] we also have 
plotted ln[s(t)] where s(t) is the mean particle size de- 
fined as 



s(t) = L(t)/N(t), 



(4) 



against ln[t], and again found a straight line with slope 
equals to 1 + 2a for all a > 0. We can thus write the 
following growth law for the mean particle size: 



s(t) ~ t 



l+2a 



(5) 



IV. ANALYTICAL MODEL 

To solve the model analytically, we use the generalized 
Smoluchowski (GS) equation 

r d d 1 f°° 

di + g^ v ( x > t )\ c ( x ' t ) = -c(x,t) J K(x,y)c{y,t)dy 
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dyK (y, x - y)c(y, t)c(x -y,t), 



(6) 



where v{x, t) is the velocity with which particles grow by 
condensation. In the absence of the second term on the 



right hand side, Eq. (|6|) reduces to the classical Smolu- 
chowski (CS) equation as given in (Eq. [1]) whose dynam- 
ics is governed by the conservation of mass law Ml. The 
GS equation does not automatically describe our model 
unless the expressions for the growth velocity v(x,t), the 
collision time r, and the kernel K(x, y) are suitably spec- 
ified. To obtain a suitable expression for v(x, t), it is 
worthwhile to recall the definition of the mean growth 
velocity which is defined as 



average growth velocity 



net growth size Ax 
elapsed time At 



(7) 



According to rule (ii), the net growth of a particle of size 
x between collisions is Ax — ax. On the other hand, 
a simple dimensional analysis of Eq. ^ reveals that 
the inverse of J K(x,y)c(y,t)dy is the collision time 
At = t(x) during which the growth ax takes place. The 
mean growth velocity between collision therefore is 



v(x, t) 



t(x) 



<>.z I dyK(x,y)c(y,t). (8) 
/o 



On the other hand, the rule (i) says that a particle col- 
lides with any particle in the system, irrespective of their 
size with an equal a priori probability. We therefore 
should choose the collision kernel independent of the size 
of the colliding particles i.e., we choose constant collision 
kernel or 



K(x,y)=2, 



(9) 



for convenience. 

To check our results from the numerical simulation, we 
now incorporate the k th moment defined as 



M k {t) 



POO 

/ x k c(x,t)dx, with k > 0, 



(10) 



together with Eqs. §E§ and © in Eq. © to write the 
rate equation for Mfc(f) in the closed form 



dM k 
dt 



E 

r=0 



M r M k -r + 2(ak-l)MoM k , (11) 



for integer k value. We can readily solve it for the first 
moment M\(t) = L(t) to give 



L(t) = (1 + JV(O)t) 



2d 



(12) 



In the long time limit it grows following the same relation 
as in Eq. (|12|) which confirms a perfect matching with 
our numerical simulation. On the other hand, solving 
Eq. pT|) with n = we find that 



N(t) 



N(0) 



i + N(o)r 



(13) 



and hence asymptotically the total number N decays fol- 
lowing the same power law 



N(t) - r 



(14) 
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with the same exponent as that of the CS equation. This 
is consistent with the assumption that the new particles 
are not nucleated. A similar temporal behavior has also 
been observed in other theories and experiments [l|, fl6| . 
One can use the two solutions for L(t) and N(t) in Eq. 
(|4j) to find the growth law for the mean particle size s(i) 



s(t) ~ (l + JV(0)t) 



l+2a 



(15) 



Therefore, one can immediately find that in the long-time 
limit the mean particle grows algebraically following the 
same relation as in Eq. ([5]) which further confirms that 
the GS equation together with the constant kernel and 
the growth velocity do describe the model in question. 



V. EXACT SOLUTION 

To solve GS equation exactly, we define the Laplace 
transformation tp(p, t) of c(x,t) by 



i/>(p,t) = / dxe~ px c(x,t), (16) 
Jo 

and its inverse transform to obtain c(x, t) by 

C(M) = 2^' dpe^(p,t), (17) 



7 — ioo 



with > 7. Differentiating Eq. (Tl7o|) with respect 

to t and using Eq. ([6]), we find that tp(p,t) obeys the 
following nonlinear partial differential equation 



dip(p, t) 



d 



0j + 2N(t)[l-ap—]^(p,t)=^ 2 (p,t), (18) 



where, of course 



N(t) = ij>(fl,t) = M (t). 



(19) 



Consequently, we need to solve Eq. (fTS]) subject to the 
initial condition 

/>oo 

lKp,0)= / dxe-^c{x,0) = f(p), (20) 
Jo 

We incorporate the solution for N(t) from Eq. (fi~3|) in 
Eq. (fl8|) and then re-write it as follows 

dxp(p, t) / 2aN{Q)p \ dip(p, t) 



dt 



l + N(0)tJ dp 

^-Trm-t^- (21) 



To solve this equation we use the method of characteristic 
in which one usually writes 



dip dip dt dtp dp 
ds dt ds dp ds 



(22) 



Comparing the above two equations we get 
dt(s) 



ds 
dp(s[ 
ds 



= 1, 



2aN(0) 
l+N{0)t 



p(s). 



(23) 



The quantity ip thus evolves following the ordinary non- 
linear equation 



2N(0) 



dM^_ _ , 2 

ds V 1 + N(0)t(s) 



4>. 



(24) 



Solving Eqs. ([23| subject to initial data 

t(s = 0) = 
p(s = 0) = po 

rP(s - 0) = f( Po ) (25) 

yields 

t = s 

Po =p{l + N(0)s) 2a . (26) 

We can transform Eq. ([M]) into a linear equation by 
setting 



to obtain 



= i 

X 



d X 2N(0) 



ds 1 + N(0)s 



X 



(27) 



(28) 



We solve it by using the integrating factor / = (1 + 
iV(0)s)~ 2 and find that 

X(s,p ) = (1 + iV(0).s)((l + N(0)s) X (0,po) - s) , (29) 
and hence using it back in Eq. (|27|) we get 
f(p(l + N(0)t) 2a ) 



(30) 



Note that, according to the definition of our model, we 
are interested in systems containing a large number of 
chemically identical particles of unit size and hence with- 
out loss of generality we may choose the mono-disperse 
initial condition 



c(x,0) = 5{x - 1), 



which gives 



f(p) = 



and AT(0) = 1. 



(31) 



(32) 



To find c(x,t) we substitute Eq. (|30|) in the definition 
of the inverse Laplace transform [(Eq. (|17p] and then a 
short calculation yields 



c(cc, t) 



t (i+ty 



(l + t)a+2a( 1+t )(l+t)-*»a 



V t > 0. (33) 



We can readily see that in the limit a — > 0j we get the 
well-known solution of the CS equation 0, \ 
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VI. SCALING DESCRIPTION 

Finding scaling or self-similar solutions is, more often 
than not, of utmost importance. These are essentially 
the solutions in the long-time limit where the distribu- 
tion function c(x,t) takes a simple universal form, in the 
sense that, it is independent of initial conditions. Most 
experimental systems do evolve to the point where such 
behavior is reached. Taking the limit t — > oo and using 
the identity 



lim 



1\» 



K) 



in Eq. (33J we get 



c(x,t) ~t-( 2+2tt W ' 



(34) 



(35) 



The structure of the solution given above is highly in- 
structive as it satisfies 



c(6 (1+2q) x, bt) = g(b)c(x,t), 



where 
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(b) = b^ 2+2a \ 



(36) 



(37) 



It implies that if we increase the units of measurement of 
a; by a factor 6( 1+2o! ) and that of time by a factor of 6, the 
numerical value of c(x, t) is decreased by a factor of g(b). 
The existence of scaling in the present case means that a 
plot of c(x, t)/t _/3 vs x/t z collapses into one graph for all 
initial conditions. Mathematically, a solution of this kind 
is called scale invariant and it implies that the system 
lacks a characteristic length scale. A further testament 
to the fact that the solution given by Eq. ([33)1 exhibits 
scaling is that it has exactly the same form as in Eq. ([2]) 
and hence comparing the two we can extract the scaling 
exponents 

(3 = 2 + 2a, z = l + 2a, (38) 
and the scaling function 

0(0 - e"*. (39) 



We can now substitute Eq. (I35[) in Eq. (jlQ|) to obtain 
the asymptotic solution for the k th moment 



where, 



M k {t)~t~^ k \ 



j(k) = 1 - (2a + l)fc. 



(40) 



(41) 



We thus see that the exponent 7(A) is linear in k which 
means that there exists a constant gap between expo- 
nents of consecutive k, and hence we can define the mean 
cluster size 



s(t) 



M k (t) 
Af k _i(t) 



for integer k > 1, (42) 



The mean particle size s(t) therefore grows algebraically 
with an exponent equal to 



7(A) - 7(* - 1) = 1 + 2a, 



(43) 



which is exactly the same as in Eq. ([5|) obtained by 
numerical simulation. Using Eq. |5J in Eq. ([2]), we 
can therefore write yet another widely used form of the 
scaling anzatz, 



z{x,t)~s{t)- e 4>(x/s{t)), 



(44) 



with the mass exponent 6 — 0/z, and hence according to 
Eq. (J38]) we get 



2 + 2a 
1 + 2a 



(45) 



The expression for the exponent j(q) reveals that the 
moment M q (t) becomes time independent if 



1 + 2a 



V a > 0, 



(46) 



since j(q) = 0. Incorporating the value of q in Eq. ([45 
we immediately find that 



1 + 9, 



(47) 



which is always less than 2 for all a > 0. We can re- 
cover the classical value = 2 of aggregation without 
condensation by setting a = 0. 



VII. DISCUSSION AND SUMMARY 

One may solve the GS equation for other growth ve- 
locities following the same method. Culille and Sire, in 
fact, solved the GS equation for v — 1 and v = x exactly 
and found 



c(x, t) 



f 2 hit 



c(x, t) 



(48) 



respectively for constant kernel K(x,y) = 1 [l5j. As 
we have the exact solutions for three different growth 
velocities, we find it worthwhile to compare them and 
see their dissimilarities. First of all, note that one can- 
not choose v = 1 and v — x and claim them to de- 
scribe the growth velocities between collisions since the 
elapsed time At during which the net growth occurs has 
not been chosen to be equal to the mean collision time 
t = 1/ J K(x,y)c(y,t)dy. Second, neither of the two 
solutions given in Eq. (|48[) can be expressed in the form 
of Eq. (|36|) or in the form of Eq. ([2]) and hence they 
violate scaling. Third, one finds that the A;th moment 
of both solutions (v = 1 and v = x) no longer exhibits 
power-law against time t which is a further testament to 
the violation of scaling since power-law relations for var- 
ious moments of the distribution function c(x, t) is the 
hall mark for the existence of scaling. Finally, when we 
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choose v = and v — ax/r, the dynamics of the systems 
are governed by some conservation laws and the scaling 
exponents are in fact fixed by these laws. However, there 
is no such conservation law in the case of v = 1 or v = x 
and since the scaling is violated and hence exponents are 
nonexistent. 

In this paper, we have presented an exactly solvable 
analytical model to study the kinetics of irreversible ag- 
gregation of particles growing by heterogeneous conden- 
sation with velocity v = ax/r. As a result of the ad- 
ditional growth by condensation, we found an algebraic 
growth law for the mean particle size s(t) ~ t 1+2a instead 
of a linear growth s(t) ~ t in the absence of growth by 
condensation. The size spectra of the aggregates formed 
by the condensation-driven aggregation is shown to ex- 
hibit universal scaling c(x,t) ~ t~ 6z (f)(x/t z ) with mass 
exponent 6 = 1 + q which is always less than its classical 



counterpart 9 = 2, and kinetic exponent z = 1 + 2a. We 
have found the exact and explicit solution for the scal- 
ing function ~ e~^. Interestingly, one obtain the 
same scaling function for aggregation without con- 
densation. The difference appears only in the z values 
of the scaling argument £ = x/t z since the two systems 
have different z values. In this sence, the form of the 
scaling function is universal in nature. We have shown 
that the transition to such scaling is accompanied by 
the emergence of a non-trivial conservation law, that is, 
M q (t) ~ constant. We believe that the present work will 
attract a renewed interest in the subject. The ideas devel- 
oped in this paper could be taken further by investigating 
the same model for generalized homogeneous aggregation 
kernel K{x,y) = {xy) x / 2 which we intend to do in our 
future endeavor. 
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